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ABSTRACT 

We present a Chandra observation of the central (r < 200 kpc) region of the cluster of galaxies 2A 0335+096, 
rich in interesting phenomena. On large scales (r > 40 kpc), the X-ray surface brightness is symmetric and 
slightly elliptical. The cluster has a cool dense core; the radial temperature gradient varies with position angle. 
The radial metallicity profile shows a pronounced central drop and an off-center peak. Similarly to many 
clusters with dense cores, 2A 0335+096 hosts a cold front at r m 40 kpc south of the center. The gas pressure 
across the front is discontinuous by a factor Ap = 1.6 ± 0.3, indicating that the cool core is moving with 
respect to the ambient gas with a Mach number M « 0.75 ± 0.2. The central dense region inside the cold 
front shows an unusual X-ray morphology, which consists of a number of X-ray blobs and/or filaments on 
scales > 3 kpc, along with two prominent X-ray cavities. The X-ray blobs are not correlated with either the 
optical line emission (Ha+[Nll]), member galaxies or radio emission. Deprojected temperature of the dense 
blobs is consistent with that of the less dense ambient gas, so these gas phases do not appear to be in thermal 
pressure equilibrium. An interesting possibility is a significant, unseen non-thermal pressure component in the 
inter-blob gas, possibly arising from the activity of the central AGN. We discuss two models for the origin of 
the gas blobs — hydrodynamic instabilities caused by the observed motion of the gas core, or "bubbling" of 
the core caused by multiple outbursts of the central AGN. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (2A 0335+096) — X-rays: 
galaxies — cooling flows 



1. INTRODUCTION 

The cores of clusters contain great dynamical complexity. 
The presence of gas cooling (Fabian 1994), relativistic plasma 
ejection from the central galaxy (McNamara et al. 2000, 
Fabian et al. 2000), the action of magnetic fields (Taylor, 
Fabian & Allen, 2002), possible thermal conduction (Voigt 
et al. 2002), local star-formation (Allen et al. 1995) and sub- 
stantial masses of cold (T < 30 K) molecular gas (Edge 200 1 ) 
all complicate the simple hydrostatic picture. 

The launch of Chandra and XMM-Newton have allowed 
the cores of clusters to be studied in unprecedented detail spa- 
tially and spectrally. These advances have led to a number 
of important results (e.g., Peterson et al. 2001; Forman et al. 
2002, McNamara 2002). The very restrictive limits on low- 
temperature gas from XMM-Newton RGS spectra are forcing 
us to re-examine some of the basics paradigms of the cooling 
flows (e.g. Kaastra et al. 2001; Peterson et al. 2001, Tamura 
et al. 2001, Molendi & Pizzolato 2001). The spatial reso- 
lution of Chandra can be used to study the brightest, nearby 
cooling flows to advance greatly our understanding of the pro- 
cesses occurring in these regions. In this paper we present the 
Chandra observation of the Zwicky cluster Zw 0335.1+0956 
whose properties may help us to address some of the above 
issues 

The Zw 0335.1+0956 was first detected as a strong X-ray 
source by Ariel-V (Cooke et al. 1978) and we therefore use 
the traditional identification of 2A 0335+096 in this paper. 
2A 0335+096 is among the brightest 25 clusters in the X- 
ray sky (Edge et al. 1990). The presence of a cooling flow 
was noted by Schwartz, Schwarz, & Tucker (1980) and its X- 
ray properties have been studied extensively over the past two 
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decades (Singh, Westergaard & Schnopper 1986, 1988; White 
et al. 1991; Sarazin O'Connell & McNamara 1992; Irwin & 
Sarazin 1995, Kikuchi et al. 1999). The optical properties of 
the central galaxy in 2A 0335+096 have been studied in de- 
tail (Romanishin & Hintzen 1988; McNamara, O'Connell & 
Bregman 1990) and the strong, extended optical line emis- 
sion marks this system out as an atypical elliptical galaxy 
but a prototypical central galaxy in a cooling flow. A deep, 
multi-frequency radio study of 2A 0335+096 (Sarazin, Baum 
& O'Dea 1995) shows a weak, flat-spectrum radio source co- 
incident with the dominant galaxy which is surrounded by 
an amorphous steep-spectrum 'mini-halo'. The tentative de- 
tection of HI absorption (McNamara, Bregman & O'Connell 
1990) and firm detection of CO emission (implying 2x109 
Mq of molecular gas) and IRAS 60 /im continuum (Edge 
2001) further highlight this cluster as one for detailed study. 
The implied mass deposition rate from the CO detection is 
low (<5 M© yr -1 ) if the cold molecular gas found is de- 
posited in the cooling flow. 

We use H = 70 km s" 1 kpc"\ Q = 0.3, and A = 0.7, 
which imply a linear scale of 0.7 kpc per arcsec at the distance 
of 2A 0335+096 (z = 0.0349). Unless specified otherwise, 
all the errors are at 90% confidence level for one interesting 
parameter. 

2. DATA PREPARATION 

2A 0335+096 was observed on 06 Sept 2000 with the Ad- 
vanced CCD Imaging Spectrometer (ACIS) using the back- 
illuminated chip S3. The total exposure time is 21.4 ksec. In 
this paper we concentrate on the bright central r < 200 kpc 
region of the cluster which lies fully within the ACIS-S3 chip. 

Hot pixels, bad columns, chip node boundaries, and events 
with grades 1, 5, and 7 were excluded from the analysis. 
We cleaned the observation of background flares following 
the prescription given in Markevitch et al. (2003). Because 
chip S3 contains the bright cluster region, we extracted a light 
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FIG. 1. — Photon image of the central region (ACIS-S3) of 2A 0335+096 in the 0.3-9 keV energy band. The image is background subtracted, vignetting 
corrected, and binned to 1" pixels. The image shows wealth of substructures in the innermost 50 kpc region. Moreover, it shows two significant surface 
brightness depressions (cavities) on opposite sides of the X-ray peak in the direction North-West South-East. We also notice a surface brightness edge to the 
South. 



curve for the other backside-illuminated chip S 1 in the 2.5- 
6 keV band, which showed that most of the observation is 
affected by a mild flare. Because the region of the cluster 
under study is very bright, we chose a less conservative flare 
cleaning than is recommended, excluding the time intervals 
with rates above a factor of 2 of the quiescent rate (taken from 
the corresponding blank-sky background dataset, Markevitch 
2001 3 ) instead of the nominal factor 1.2. This resulted in the 
final exposure of 13.8 ks. During the accepted exposure, the 
frontside-illuminated chips did not exhibit any background 
rate increases, therefore the residual background flare is of the 
"soft" species and can be approximately modeled as described 
by Markevitch et al. (2003). We fitted the excess background 
in chip S 1 above 2.5 keV with the flare model and included it, 
rescaled by the respective solid angle, to the cluster fits. The 
main component of the background was modeled using the 
above-mentioned quiescent blank-sky dataset, normalized by 
the rate in the 10-12 keV band (which can be done because 
the residual flare in our observation is "soft" and does not 
affect the high-energy rate). The addition of the flare compo- 
nent has a very small effect on the results, because the cluster 
core is so bright. We therefore ignored the possible (not yet 
adequately studied) spatial nonuniformity of the flare compo- 

3 http://asc.harvard.edu/ "Instruments and Calibration", "ACIS", "AOS 
Background" 



nent. 

3. X-RAY IMAGING ANALYSIS 

Fig.[2shows a background subtracted, vignetting corrected 
image of the central region of the cluster, extracted in the 0.3- 
9 keV energy band and binned to 1" pixels. The cluster X-ray 
surface brightness appears to be regular and slightly elliptical 
at r > 70" (50 kpc). To the South, at about 60" (42 kpc) 
from the X-ray peak, we notice the presence of a sharp surface 
brightness edge similar, to those observed in other clusters of 
galaxies (e.g., Markevitch et al. 2000, Vikhlinin, Markevitch, 
& Murray 2002a; Mazzotta et al. 2001). Finally, the X-ray im- 
age shows complex structure in the innermost rj 50" (35 kpc) 
region. Below we study these features in detail. 

3.1. Core structure on larger scale 

To investigate the cluster spatial structure, and to compare 
it with previous studies, we start with the classic approach of 
fitting the surface brightness with a /3-model (Cavaliere, & 
Fusco-Femiano, 1976). Consistent with previous results, we 
find that this model does not provide an acceptable fit. Hence, 
we use a double 2-dimensional /3-model defined by: 

^ f narrow "F /extended- (1) 

The narrow and extended components of eq.[^are given by: 



f(r,r c ,f3,x ,y ,e,6) = A x l+(— ) 



-3/3+1/2 



(2) 
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Table 1 

Result of fitting the central 250" cluster X-ray image with a double 2D-elliptical-/3-model 



Component 


x, y 


A r c 





e 




(J2000) 


(cnt arcsec -2 ) (arcsec) 




(deg) 


Narrow 


(03:38:40.9; +09:57:57.5)1 


13.2 ± 0.2 85 ±1 


2.45 ±0.04 0.12 ±0.008 


33 ±6 


Extended 


(03:38:40.4; +09:58:31.7)1 


5.05 ±0.06 64 ±1 


0.571 ±0.004 0.13 ±0.01 


34 ± 4 



fThe error of the central position is 2' 




Right Ascension 



FIG. 2. — Comparison between the X-ray image of 2A 0335+096 and 
the best fit double /3-model. Continuous lines: isointensity contour levels of 
the cluster image shown in Fig.fTlafter adaptive smoothing. The image was 
smoothed using a Gaussian kernel with a S/N ratio of the signal under the 
kernel > 3. The levels are spaced by a factor \f2. Dashed lines: isointensity 
contour levels of the best fit double /3-model. The levels are the same as for 
the X-ray image. 
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FIG. 3. — Comparison between the surface brightness radial profiles ex- 
tracted from the X-ray image and the best fit double /3-model. To extract the 
profiles we used concentric ellipses whose center, ellipticity and the orienta- 
tion where equal to the best fit values of the narrow component of the double 
/3-model. Points with errorbars and lines correspond to the profiles extracted 
from the X-ray image and the double /3-model, respectively. 



where 

r(x,y) = V ^^ 1 ^' + l £ ; (3) 

and 

Xnew = (x- x ) cos(0) + (y - y D ) sin(6>). (4) 

To perform the fit we use the SHERPA fitting program 4 . 
The cluster region with r > 250" and all the stro ng p oint 
sources, together with the regions defined later in fl3.3l are 
masked out. Because the image contains a large number of 
pixels with few number counts (< 20 cnt) we use the Cash 
(1979) statistics which, unlike the \2 one, can be applied re- 
gardless of the number of counts per bin. The best fit values 
with their 90% errors are reported in Table 1 . 

Unlike the \ 2 statistics, the magnitude of the Cash statistics 
depends upon the number of bins included in the fit and the 
values of the data themselves. This means that there is not 
an easy way to assess from the Cash magnitude how well the 
model fits the data. To have a qualitative idea of the goodness 

4 http://asc.harvard.edu/ "Data Analysis", "Introduction" 



of the fit, we generate a mock X-ray image using the best fit 
values of the double /3-model and we apply the Poisson scatter 
to it. Consistent with previous studies, we find that globally 
the best fit double /3-model describes very well the overall 
cluster structure. This is evident in Fig.|2]where we show the 
isocontour levels of the model image overlaid to the same iso- 
contour levels of the cluster image after adaptive smoothing. 
To better illustrate the quality of the fit we also extract the X- 
ray surface brightness radial profiles in concentric elliptical 
annuli from both the cluster and the mock images as shown in 
Fig- El The center, the ellipticity and the orientation of each 
ellipse are set to the best fit values of the narrow component 
of the double /3-model. The figure clearly confirms the good 
quality of the fit. 

The results of the fit reported in Table 1 show that the two 
spatial components have very different slopes, in particular 
the narrow component is much steeper than the extended one. 
We also notice that the profile of the narrow component de- 
rived from the Chandra data is significantly steeper that the 
one obtained from the ROSAT PSPC (Kikuchi et al. 1999). 
The large value associated with the narrow component in- 
dicates that the cluster contains a very compact and dense gas 
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FIG. 4. — Upper panel) Cluster surface brightness profile extracted from 
the 150° to 210° sectors. The continuous line represents the best-fit of the 
discontinuous two-power law density model. Lower panel) Deprojected tem- 
perature profile from the 150° to 210° sector. The continuous line is the 
best-fit of the single power law model. 



Table 2 

x-ray front density, temperature, and pressure best 
fit values "("(errors 90% c.l.) 



Quantity (Q) £1 £2 A r,„ mp (arcsec) 

Density (n) -1.36 ± 0.04 -1.20 ± 0.02 1.6 ± 0.02 56.7 ±1 
Temperature (T) +0.79 ± 0.03 +0.67 ±0.20 1.0 + 0.20 
Pressure (P) -0.57 + 0.05 -0.53 + 0.20 1.6 + 0.30 
fEach quantity is modeled using two power-laws with a jump: 

q I (^/^jump*) 1 T *Tjump 

\ Aq\T I T jump*) 2 T ^ l"jump 



Table 3 

Spatial and spectral properties of the X-ray blobs 



Region 


n 


ri 


net cnts 


n e 


^[0.6-9] 


tc 




(arcsec) ( 


arcsec) ([0.6-9] keV) (10 


-2 

'cm 


' 3 ) 10 41 erg s- 1 


(108 5 


1 


9.0 


4.8 


796 


8.9 


5.8 


2.0 


2 


7.5 


5.6 


1092 


10.0 


8.2 


1.9 


3 


5.8 


4.6 


724 


12.7 


6.4 


2.1 


4 


6.6 


3.7 


582 


9.4 


4.6 


2.4 


5 


7.4 


5.1 


405 


6.4 


2.7 


2.4 


6 


8.7 


5.1 


684 


7.6 


7.5 


4.1 


7 


4.4 


4.3 


223 


8.3 


1.8 


3.3 


8 


5.2 


5.2 


312 


7.2 


2.4 


2.8 



component with a quite sharp boundary (this is also clearly 
visible in Fig. [3. On the other hand, the slope of the ex- 
tended component derived from the Chandra data is signif- 
icantly lower that the value derived from ROSAT (see e.g. 
Kikuchi et al. 1999, Vikhlinin, Forman, & Jones 1999). Be- 
cause the Chandra result is obtained fitting the cluster sur- 
face brightness in a much smaller radial range (r < 250") 
than the ROSAT one (r < 1200"), this indicates that 2A 
0335+096 contains one or more radial breaks at r > 250", 
similar to the ones observed in other clusters like e.g. Hydra- 
A (David et al. 2001) and Perseus (Fabian et al. 2000). 

Finally we note that, although the ellipticity and the orienta- 
tion are similar, the centroid of the two spatial components are 
offset by r* 35" with the narrow component centroid being 
further south than the extended one. 

Although spatial structure of the cluster is globally well de- 
scribed by a simple double /3-model, there are two regions 
whose structure departs significantly from the best fit model 
above: i) a 60° wide sector to the south; ii) the innermost 
r < 50" cluster region. Below we study in details the spatial 
structure of these regions. 

3.2. X-ray "front" 

A closer look at Fig. [T] indicates that 2A 0335+096 hosts 
a surface brightness drop at r rj 50 — 60" from the X-ray 
peak to the south. This X-ray drop spans a sector from 150 to 
210 deg and most closely resembles similar features in other 
clusters with dense cores and regular morphology on larger 
linear scales, such as RXJ1720+26 and A1795 (Mazzotta et 
al. 2001; Markevitch, Vikhlinin, & Mazzotta 2001). To better 
visualize this feature and to study its nature, we extract the 
cluster X-ray surface brightness profile within the above sec- 
tor, using elliptical annuli concentric to the brightness edge. 
The center, the ellipticity and the orientation of each ellipses 
are set to the best fit values of the narrow component of the 



double /3-model. The profile is shown in the upper panel of 
Fig. @] As extensively shown in previous analysis of cold 
fronts in clusters of galaxies, the particular shape of the ob- 
served surface brightness profile indicates that the density dis- 
tribution is discontinuous at some radius r Jump . 

To estimate the amplitude and the position of the density 
jump, we fit the profile using a gas density model consisting 
of two power laws with index ol\ and a 2 and a jump by a 
factor A n at the radius rj ump : 

Q = \ ( r I > 'jump) 1 r < r jump ; ^\ 

|^ A n (r I v j urn p) 2 v > rj U mp • 

The fit is restricted to radii r = 30 — 200" to exclude the 
central region with irregular structure. We assume spherical 
symmetry. The best fit values together with the 90% confi- 
dence level errors are shown in Table 2. The best fit model 
(which provides a good fit to the data) is shown with a dashed 
line in the upper panel of Fig.|4] 

3.3. Core structure on smaller scales 

To study the cluster morphology on smaller scales we adap- 
tively smooth the cluster image as shown in Fig. [5] The image 
was smoothed using the CIAO 5 tool CSMOOTH with a Gaus- 
sian kernel and a S/N ratio of the signal under the kernel > 3. 

The image clearly shows a very complex structure with a 
number of blobs of X-ray excess emission. The unsmoothed 
image shows that these blobs are not necessarily distinct but 
may form a filamentary structure. Nevertheless, for the pur- 
pose of our analysis we assume that the X-ray blobs are el- 
lipsoid. We select 8 regions as shown in Fig. [5] Each ellipse 
indicates the region where the surface brightness drops by al- 
most a factor two with respect to the peak of its respective 
X-ray blob. To measure the volume occupied by each blob 

5 http://asc.harvard.edu/ "Data Analysis" 
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Table 4 

Comparison of Spectral Fits for the innermost 250'' 



REGION 



Model 


N H (10 21 cm" 2 ) 


Ti (keV) 


T 2 (keV) 


z/z 


M (M Q yr- 2 ) 


X 2 /d.o./. 


TBABS x MEKAL 


1.6 

(i.6±8:S) 


(2.6+r |) 




0.59 
(0-54±g;g*) 




843/352 
(484/297) 


TBABS x [MEKAL+MEKAL] 


-i 71 +0.06 
' 1 -0.05 

(1.7ltg : g?) 


o O+0.9 
O.o_ 2 

C3-0±g;|) 


1.4±g;S 

(i.4±S:i) 


Co.53±g : gS) 




643/350 
(393/295) 


TBABS x [MEKAL+MKCFLOW] 


2.0±° ;1 

(2.olS:i) 


2.90+°l 
(2.78lg;° 6 6 ) 




0.68lg : g* 

(o.6olg:°*) 


(192111) 


658/351 
(387/296) 


TBABS x [MEKAL 
+ZTB AB S (MKCFLOW) ] 


o.9±g;lt 

(i.3tS:l)t 


(2.8l° ;D 




0.68l° : g 


266+^ 
(249lg) 


658/351 
(393/296) 



tlntrinsic absorption 

Values in parentheses are from fits with the excluded 1.4-2.2 keV energy interval. 




FIG. 5. — Zoom in of Fig.fTlafter adaptive smoothing made using a Gaussian 
kernel with a S/N ratio of the signal under the kernel > 3. The green ellipses 
labeled from 1 to 8 identify 8 X-ray bright extended regions or blobs. The 
two blue ellipses, labeled 9 and 10, identify the Western and the Eastern X- 
ray holes or cavities. The white square and cross, indicate the position of 
the cluster D galaxy and its companion, respectively. The wh ite X symbol 
indicate the position of the other cluster central galaxy see (Fig. II II . 



we make the simple assumption that the radius of the ellip- 
soid along the line of sight is equal to the minor radius of the 
ellipse r 2 . In Table 3 we report the major and the minor radii 
corresponding to each region shown in Fig. [5] The dimen- 
sions of the X-ray blobs span from 4" to 9" (from 2.8 to 
6.3 kpc). For each blob we calculate the mean electron den- 
sity using the following simple procedure: using the method 
described later in § |4j we extract the spectra from both the 
blob region and from a cluster region just outside it. Using 
the latter spectrum as background, we fit the blob spectrum 
with a thermal model, hence deriving the emission measure. 
Finally, the gas density is calculated from the emission mea- 
sure under the simple assumption that it is constant within the 
blob. The net number of counts and the mean gas densities 
obtained using the above procedure are reported Table 3. For 



completeness in the same table we add the unabsorbed lumi- 
nosity in the (0.6- 9) ke V band in the source rest frame and the 
cooling time (see tj4.5l below). 

Beside the X-ray blobs, we also note a prominent X-ray 
hole about 30" (21 kpc) North- West from the X-ray peak, and 
a less prominent one at « 40" (38 kpc) to the East. To check 
the statistical significance of the Eastern hole, we extract the 
count rates from a circular region in the hole (r = 6.5"; see 
Fig. |5} and the count rates in two similar regions to the East 
and to the West of the depression along the isocontour level 
defined by the best-fit double /3-model. The absorbed flux 
from the hole is « 20% lower than the flux from the neigh- 
boring regions at a 4.5er confidence level. As shown later in 
||5] the cluster image may also present other X-ray holes at 
larger radii. The previous identification of the Eastern hole 
from ROSAT HRI observations with absorption by a cluster 
member galaxy by Sarazinetal. (1992, 1995) is not supported 
with these higher resolution images (see Fig.|5Jl. 

4. SPECTRAL ANALYSIS 

We use the CIAO software to extract spectra. The spectra 
are grouped to have a minimum of 20 net counts per bin and 
fitted using the XSPEC package (Arnaud 1996). The position- 
dependent RMFs and ARFs are computed and weighted by 
the X-ray brightness over the corresponding image region us- 
ing "calcrmf" and "calcarf" tools 6 . To account for the ab- 
sorption caused by molecular contamination of the ACIS op- 
tical blocking filters 7 , all the spectral model we use are mul- 
tiplied by the ACISABS absorption model 8 that models the 
above absorption as a function of time since launch. The pa- 
rameter TDAYS is set to 410, while the others are left to the 
default values, namely NORM= 7.22 10~ 3 , TAUINF= 0.582, 
tefold= 620, nC= 10, nH= 20, nO= 2, nN= 1. 

Spectra are extracted in the 0.6-9.0 keV band in PI channels. 
At present, the ACIS response is poorly calibrated around the 
mirror Ir edge, which results in frequently observed signifi- 

6 A. Vikhlinin 2000, i http://asc.harvard.edu/ "Software Exchange", "Con- 
tributed Software"). 

7 http://asc.harvard.edu/ "Instrument & Calibration", "ACIS", "ACIS QE 
Degradation" 

8 Chartas & Getman 2002, http://www.astro.psu.edu/users/chartas/xcontdir/xcont.html 



6 



MAZZOTTA ET AL. 




FIG. 6. — Temperature map of 2A 0335+096 overlaid on the 0.3-9 keV ACIS brightness contours. Left panel) This temperature map was obtained by smoothing 
the images using a Gaussian filter with cr = 15". The green ellipse corresponds to temperature map on the right panel. Blue holes indicate the position of the 
point sources that were masked out. The isocontours are obtained from the image after a Gaussian smoothing with a = 15". Levels are spaced by a factor 
\/2. Right panel) Zoom in, high resolution temperature map from the cluster inside the green ellipse on the left panel. This temperature map was obtained by 
smoothing the images using a Gaussian filter with a = 4". Isolevels are the same as in the continuous lines in Fig|2| The black ellipses labeled from 1 to 10 
identify the regions defined in Fig|5] The blue square, cross, and X symbols indicate the position of the galaxies as in Fig. 151 



cant residuals in the 1.4-2.2 keV energy interval and high \ 2 
values. We tried excluding this energy interval and found that 
the best-fit parameter values and confidence intervals do not 
change considerably while the x 2 values become acceptable, 
as will be seen below. We therefore elected to use the full 
spectral range. 

4.1. Average Temperature 

To determine the average spectral properties of this obser- 
vation, we extract the overall spectrum from the innermost 
circular region of r = 250" (175 kpc) centered on the X-ray 
peak and fit it with different models as listed in Table 4. In 
the following the absorption of X-rays due to the interstellar 
medium has been parametrized using the Tubingen-Boulder 
model (TBABS in XSPEC v. 11.1; Wilms, Allen & Mc-Cray 
2000). Furthermore the metallicity refers to the solar photo- 
spheric values in Anders & Grevesse (1989). As a first at- 
tempt, we use an absorbed single temperature thermal model 
(TBABS XMEKAL; see e.g. Kaastra, 1992; Liedahl, Osterheld, 
& Goldstein, 1995; and references therein). We find that this 
model does not reproduce the spectrum. We also use an ab- 
sorbed two-temperature model (TBABS x [MEKAL+MEKAL]) 
with the metallicity of the two components linked together. 
The hydrogen equivalent column density and the temper- 
ature, the metallicity, and the normalization of each ther- 
mal component are left free to vary. This model provides 
a much better fit. In particular, if we exclude the energy 
band near the Ir edge from the spectral analysis, we find 
that the model is statistically acceptable (see Table 4). We 
find that the equivalent column density is consistent with 
the galactic value (N H = 1.78 x 10 21 cm~ 2 ; Dickey & 
Lockman 1990). The temperatures of both the hot and 
the cool components, as well as the metallicity, are con- 
sistent with the values obtained from the combined analy- 
sis of ASCA GIS and SIS spectra extracted from the inner- 
most r = 120" region (Kikuchy et al. 1999). We finally 



fit the spectrum using two models with cooling flows. The 
first is an absorbed thermal model plus a cooling flow com- 
ponent (tbabs x[mekal+mkcflow]). In the latter we as- 
sumed that the cooling flow component itself is intrinsically 
absorbed by uniformly distributed amount of hydrogen at the 
cluster redshift (tbabs x [mekal+ztbabs(mkcflow)], in 
XSPEC v. 11.1). The higher temperature and metallicity of 
the cooling flow component are fixed to be equal to the ones 
of the thermal component. In the first cooling flow we leave 
the absorption free to vary. In the intrinsically absorbed cool- 
ing flow model we fix the overall absorption to the Galactic 
value and let the intrinsic Nu free to vary. We find that both 
models give temperatures, metallicities and the mass deposi- 
tion rates consistent within the errors. It is worth noting that 
the cooling flow models give a minimum \ 2 value similar to 
the two-temperature model which, as said before, is quite high 
compared to the degree of freedom (d.o.f). 

4.2. Temperature Map 

To determine the temperature map we used the same tech- 
nique used for A3667 (Vikhlinin, et al. 2001a; Mazzotta, 
Fusco-Femiano, Vikhlinin 2002). We extracted images in 10 
energy bands, 0.80 - 1.07 - 1.43 - 1.68 - 1.82 - 2.16 - 
2.35 - 2.88 - 3.75 - 5.83 - 9.00 keV, subtracted background, 
masked out point sources, divided by the vignetting maps and 
smoothed the final images. We determined the temperature by 
fitting the 10 energy points to a single-temperature absorbed 
MEKAL plasma model with metal element abundance fixed 
at Z = 0.6 Zq. We left as free parameters the equivalent 
absorption column density, the temperature, and the normal- 
ization. 

To study the cluster temperature structure at larger radii 
(where the cluster is less bright), we derive first a relatively 
low spatial resolution temperature map. This is obtained by 
smoothing the 10 images above using a Gaussian filter with 
a = 15". The "low-resolution" temperate map, overlaid on 
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FIG. 7. — Emission weighted temperature, metallicity, and equivalent col- 
umn density radial profiles of 2A 0335+096 . The spectra are extracted in 
concentric elliptical annuli. Center, ellipticity and orientation of each ellipses 
is fixed to the best fit values of the narrow component of the double /3-model 
as shown in Fig.|3] The dotted crosses in the upper panel indicate the results 
obtained by deprojecting the temperature profiles in the assumption of spher- 
ical symmetry. The dotted line in the bottom panel is Njj value obtained by 
fitting a constant to the Njj profile. Error bars are la confidence level. 



the 0.3-9 kev ACIS brightness contours, is shown in the left 
panel of Fig. [6] Like the surface brightness, the projected 
temperature map is not azimuthally symmetric. The gas is 
substantially cooler in the cluster center and shows a strong 
positive temperature gradient which varies with position an- 
gle. The temperature gradient seems to be higher and lower to 
the Southern and Northern sectors, respectively. It is interest- 
ing to note that the gradient seems to be maximum somewhere 
in the south sector, where we observe the surface brightness 
edge. 

To study the the rma l properties of the X-ray blobs and cav- 
ities discussed in i|3.3l we produce a second temperature map 
of the innermost central region using a much higher spatial 
resolution. This region is identified by the green ellipse in left 
panel of Fig. [6] The "high-resolution" temperature map, over- 
laid on the 0.3-9 kev ACIS brightness contours, is shown in 
the right panel of Fig. [6] The map is obtained by smoothing 
the 10 images above using a very narrow Gaussian filter with 
a = 4". For clarity we plot ellipses numbered from 1 to 10 
corr espo nding to the 8 X-ray blobs and 2 cavities defined in 
the A3. 31 above (see also Fig. [5}- The high-resolution tem- 
perature map shows a very complex structure. In particular 
the cavities and the X-ray blobs seem to have qui te di fferent 
projected temperatures. We discuss this aspect in i)4.5l below. 

4.3. Radial temperature gradients 

We first measure the overall spectral radial properties ex- 
tracting spectra in concentric elliptical annuli. Center, elliptic- 
ity and orientation of each ellipse are fixed to the best fit val- 
ues of the narrow component of the double /3-model as shown 
in Fig. [3] The radii are choose in order to have « 20, 000 
net counts per spatial bin. Each spectrum is fitted with an ab- 
sorbed one-temperature MEKAL model. Equivalent hydrogen 
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FIG. 8. — Comparison of the temperature profiles extracted from the North- 
ern half (filled triangles) and Southern half (open circles) of the cluster. 



column density, temperature, metallicity, and normalization 
are left free to vary. We get an acceptable fit for all the annuli 
but the innermost one. The resulting temperature, metallicity, 
and Njj profiles, with error bars at la level, are shown in the 
upper, middle, and lower panel of Fig. respectively (errors 
bars are la level). The temperature profile shows a positive 
gradient within the innermost 50 — 60hjQ kpc with the tem- 
perature that goes from rj 1.5 keV to « 3.5 kev. At larger 
radii the temperature profile shows a much flatter profile con- 
sistent with being constant. Only the last bin at 150h^ kpc 
indicates a possible temperature increment at larger radii. Al- 
though the cluster is clearly not spherical symmetric we try 
to deproject the overall temperature profiles assuming it actu- 
ally is. The deprojected profile is reported as dotted crosses 
in the upper panel of Fig. [7] We notice that the deprojected 
temperature of the innermost bin is T = 1.36 ± 0.02 keV. 

The projected metallicity is consistent with the constant 
value Z w 0.5Z Q at radii > 50 — 60h^ kpc. At lower 
radii it shows strong positive gradient with a maximum of 
~ O.&Zq at AOh^Q kpc. Similar metallicity gradients have 
been observed in other clusters, e.g. M87 (Molendi & 
Gastaldello, 2001), Persus (Schmidt, Fabian and Sanders, 
2002), Abell 2052 (Blanton, Sarazin, McNamara, 2002), 
A2199 (Johnstone et al. 2002) and Centaurus (Sanders & 
Fabian 2002). For some of them, however, it has been shown 
that the gradient is an artifact induced by the inadequate mod- 
eling of the spectra with a single temperature component 
where a multi-temperature one is actually required. This ef- 
fect is known as metallicity bias (see Buote 2000). It is inter- 
esting to note, however, that if we deprojected the metallicity 
profile assuming spherical symmetry we find that the profile 
is consistent with the projected one. This either indicates that 
the spherical symmetry assumption is inadequate or that the 
gradient is real. Recently Morris & Fabian (2003) suggested 
that apparent off-center peaked metallicity profiles may form 
as result of the thermal evolution of intracluster medium with 
a non-uniform metal distribution. 
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Finally, as shown in Fig. [7] panel c), we find that the equiva- 
lent column density profile is consistent with a constant value 
of 1.98 x 10 21 cm -2 which seems to be 10% higher than the 
galactic value. 

As discussed in £14.21 the clusters temperature distribution is 
strongly asymmetric. To verify the statistical significance of 
this asymmetry, we extract the spectra using the photons from 
the Northern or the Southern sectors. Spectra are extracted 
in the elliptical annuli, with the same binning as before. We 
find that, although both the abundance and metallicity profiles 
from the North and the South sectors are consistent within 
the errors, the temperatures profiles are significantly differ- 
ent. Fig. [8] clearly shows that the temperature gradient in the 
North sector is shallower than the gradient in the South sec- 
tor. Moreover, the region of constant temperature in the South 
sector is a factor 1.3 hotter than that in the North sector. 

4.4. Temperature profile of the "front" 

In £13.21 we show that the cluster has a surface brightness 
edge similar to the cold fronts observed in other clusters of 
galaxies. One peculiar characteristic of these cold fronts is 
that they have a temperature discontinuity (jump) over the 
front itself. In this paragraph we study the thermal proper- 
ties of the X-ray "front" in 2A 0335+096. To do this we es- 
timate the 3D temperature profile of the front by deprojecting 
the temperature profile from the cluster sector whose angles 
are from 150° to 210°. The 3D temperature profile is shown 
in the lower panel of Fig.|4] We notice that the temperature 
profile has a quite steep gradient. Moreover, unlike the cold 
fronts observed in other clusters of galaxies, it does not show 
a clear temperature discontinuity at the front. At first we fit 
the temperature profile in the range from 30" to 150" with 
a single power-law model T oc 9 s . The model provides a 
good fit with S = 0.72 ± 0.02, as shown in Fig.g] We no- 
tice that the data may also be consistent with a small tem- 
perature jump at the front. To constrain the amplitude of this 
possible jump, we fit the data using a two-power law model 
with a jump At at the position of the front (see eq |3J and 
we report the result in Table 2. From the table we see that 
At = 1 ± 0.2. Using the pressure and the density models 
above, we calculate the pressure profile near the front region. 
The best fit model together with the 90% errors are reported 
in Table 2. The observed density and temperature discontinu- 
ity at the front result in corresponding pressure discontinuity 
of a factor A P = 1.6 ± 0.3. 

4.5. Spectral properties of the X-ray blobs and cavities 

As discussed in £14.21 the innermost region of 2A 
0335+096 has a very complex temperature structure. In par- 
ticular, the temperature map shows that the cavities and the 
X-ray blobs have quite different projected temperatures (see 
Fig-H). To better study the thermal properties of these struc- 
tures, we extract spectra for all the elliptical regions shown in 
the Figure and fit them with an absorbed thermal model. We 
verify if there is any excess absorption by fixing the metal- 
licity to Z=0.6 solar and letting both Nh and the temperature 
to vary freely. As we find that the resulting Nh values are 
all consistent with the cluster mean value, we fix Nh to this 
value and re-fit the model to measure the temperature. The 
projected temperatures together with the 90% error bars, are 
shown in Fig. [5] Consistently with the right panel of Fig. [5] 
we find that the central X-ray structures have different pro- 
jected temperatures. In particular we notice that region 1 and 
2 are significantly cooler than region 6 and 7. Furthermore, 
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FIG. 9. — Projected temperature from the 8 regions defined in Fig. [5] The 
two filled triangles give the projected temperature from the Western and East- 
ern holes, respectively. Error bars are 90% confidence level. 
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FIG. 10. — Upper panel) Deprojected temperature associated with each 
of the X-ray blobs region defined in Fig. [3] The dashed lines are the 90% 
eiTor measurement on the deprojected temperature of the ambient gas. Middle 
panel) Thermal pressure associated with the same blobs as in the Upper panel. 
The dashed lines are the 90% error measurement of the thermal pressure of 
the ambient gas. Lower panel) Cooling time associated with the same blobs 
as in the Upper panel. 



we notice that the holes are significantly hotter than any other 
X-ray blobs. It is worth noting that the observed higher tem- 
perature from the holes is likely to be the result of a projection 
effect. 

Projection could also be responsible for the observed dif- 
ference in the projected temperature of the X-ray blobs. To 
verify this hypothesis one should find a way to disentangle 
the contribution to the spectrum due to the blob itself from 
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that due to the gas along the line of sight. To perform this 
task we assume that the contribution to the spectra due to the 
gas along the line of sight is almost the same within the in- 
nermost 25" region (where all the blobs are). This is justified 
by the fact that the surface brightness profile within this re- 
gion does not vary more than a factor » 1.5. We then extract 
a cluster spectrum using a region just outside the blobs that 
will be used as background for the spectrum of each blob. 
The first important result of this exercise is that the final spec- 
tra clearly shows the presence of Fe-L emission line complex. 
This proves, without doubt, that the emission from these blobs 
is indeed mostly, if not all, thermal. We fit these spectra using 
a single-temperature MEKAL model leaving as a free param- 
eter only the temperature and the normalization. The result 
is shown in the upper panel of Fig. ^3 I n the same figure 
we add, as dashed lines, the 90% error temperture estimate 
of the ambient gas. This estimate cames from the fit of the 
blob-free cluster spectrum, used before as background, with 
a two-temperature MEKAL model with the metallicity of the 
two thermal components linked together. To compensate our 
ignorance of the actual 3D cluster structure, we leave the nor- 
malization of both thermal components free to vary. We as- 
sume that the temperature of the ambient gas coincides with 
the one of lower temperature component. We believe that this 
represents a good assumption as the value we find is consis- 
tent with the temperature of the innermost bin obtained using 
a spherical symmetric deprojection model (see upper panel of 
Fig.Q. 

Unlike the emission weighted temperatures, the deprojected 
ones have a much lower scatter. Furthermore, all regions, but 
the 6th, have temperatures consistent with the 3D temperature 
of the ambient gas. 

Using the blob density measurements reported in Table 3, 
it is now possible to estimate the thermal pressure associated 
with each blob. This is shown in the central panel of Fig. 1101 
In the same figure we add, as dashed lines, the 90% error 
ambient pressure. The blobs, except region 5, have thermal 
pressure values significantly higher than the ambient thermal 
pressure. This unexpected finding will be discussed below. 

Using the temperature and the density of each blob we also 
estimate their cooling times r c = 5P/2e, where P and e are 
the blob's pressure and the emissivity, respectively. These, to- 
gether with the 90% error bar, are shown in the lower panel 
of Fig. [To] We notice that all blobs have similar cooling times 
which are consistent within the errors (but see region 6). For 
practical reasons we report also the blob cooling times in col- 
umn 7 of Table 3. 

For completeness we estimate the X-ray blob unabsorbed 
luminosities in the (0.6-9.0) keV in the source rest frame, as 
shown in in column 6 of Table 3. 

5. PROPERTIES AT OTHER WAVELENGTHS 

We conclude the analysis of 2A 0335+096 by comparing 
the cluster properties in the X-rays to those at other wave- 
lengths. 

In Fig.^2 we show an HST image of the central region of 
2A 0335+096 superposed on contours of the X-ray image af- 
ter adaptive smoothing. The optical image shows one central 
D galaxy (PGC 013424, Paturel et al. 1989) with a compan- 
ion smaller galaxy 7" to the North- West and another elliptical 
galaxy at 41" t o the East. In Fig.[^the X-ray blobs, defined 
previously in fl3.3l (labeled 1 to 8 in Fig. |5J, can be easily 
identified matching the small quasi-elliptical isocontour re- 
gions. Apart from one X-ray blob coincident with the cluster 




FIG. 1 1 . — HST image of the central region of 2A 0335+096 superposed on 
contours of the X-ray image after adaptive smoothing. The circles in square 
symbols indicate the best fit central position of each of the components of 
the double /3-model. The Northern and the Southern symbols indicate the 
position of the Extended and Narrow component, respectively. The thick 
ellipses indicate the position of the Western and Eastern holes, respectively. 
The X-ray excess regions labeled 1 to 8 in Fig. |5] can be easily identified 
by the small quasi-elliptical isocontour regions. Notice that none of these 
regions have an optical counterpart. 



D galaxy none of the other blobs have optical counterparts. 
This strongly suggests that the X-ray blobs are produced by 
some kind of hydrodynamic process. 

For completeness in Fig.^2 we report the centroid position 
of the narrow and extended /3-model components. It is worth 
noting that the D galaxy position does not coincide with either 
of these centroids but rather lies in between. 

The central galaxy of 2A 0335+096 has been found to have 
a number of peculiar properties over the past 2 decades. Ini- 
tially it was classified as a Seyfert 1.5 or Liner by Huchra, 
Wyatt & Davis (1982) but more detailed optical observations 
by Romanishin & Hintzen (1988) found a very extended opti- 
cal emission line region (Ha+[Nll]) to the North-East of the 
galaxy (extended over 30 kpc) and the central regions of the 
galaxy are anomalously blue indicating recent star-formation 
(see Fig. 1121. Recently a number of authors have shown that 
some of the extended Ha features observed in the central 
region of a number of clusters correlate with similar struc- 
tures of X-ray excess (e.g. A1795, Fabian et al. 2001; RX 
J0820.9+0752, Bayer-Kim et al. 2002; Virgo, Young, Wilson, 
& Mundell 2002). 

To study the correlation between the X-ray blobs and the 
line emission in Fig.^]we show the narrow band Ha+[Nll] 
image from Romanishin & Hintzen (1988) with the X-ray iso- 
contours superimposed. The accuracy of the alignment is bet- 
ter than 3". There is no direct correlation between the X-ray 
and Ha+[Nll] structures but there is a tendency for the two 
to lie close to each other. However, the statistical significance 
of this tendency is difficult to judge. Note the X-ray blobs 
labeled 1 and 4 lie too close to the saturated star to the South 
East of the Ho+[Nll] image to detect the fainter Ha+[Nll] 
emission. 

The radio properties of 2A 0335+096 are also of note. 
Sarazin et al. (1995) find an amorphous, steep spectrum 
(a ~ -1.4) 'mini-halo' around PGC 013424 which has 
relatively low radio power compared to other similar sys- 
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FIG. 12. — Ha+[NII] image of the central 60" X 70" region of 2A 0335+096 (from Romanishin & Hintzen 1988, Fig. 3) with X-ray contours overlaid. The 
X-ray contours are the same as in Fig. II II To better visualize the faintest line features we show the Ha+[NlI] image in both high and low contrast (left and right 
panels, respectively). The square and the cross symbols indicate the position of the central D galaxy and its companion, respectively. The X-ray excess region 
labeled 1 to 8 in Fig.|5]can be easily identified by the small quasi-elliptical isocontour regions. The Western X-ray cavity is located in the upper right hand corner. 
Note that the large, bright artifact in the lower left of the image is a saturated star. 




FIG. 13. — Percentile variation of the cluster X-ray image (after adap- 
tive smoothing) with respect to the best fit double /3-model superposed on 
contours of the 1.515 GHz radio image (from Sarazin, Baum, & O'Dea, 
1995). The darker and lighter regions indicate negative and positive variation 
> |20|%, respectively. The arrows indicate the position of the Western and 
Eastern holes and the position of the Front. The square and the cross symbols 
indicate the position of the central D galaxy and its companion while the X 
sym bol indicates the position of the third cluster galaxy in the HST field (see 
Fig. 1111 . Finally the circles in squares symbols indicate the best fit central 
position of each of the components of the double /3-model. 



tems. Fig. ^] shows the image of the percentile variation 
((S c ; — S mo£ j)/S c ;) of the cluster X-ray surface brightness 
S c ; with respect to the best fit double /3-model £ mo d- The 
isocontours of the VLA 1.5 GHz C+D-array image (resolu- 
tion 16") from Sarazin et al. (1995) are superimposed. The 
image clearly shows the position of the X-ray front together 
with the X-ray blob s an d the Western and Eastern holes ear- 
lier studied in A3. 21 and fl3.3l respectively. We notice that the 
mini-radio halo is well within the radius defined by the posi- 
tion of the X-ray front. Unfortunately the spatial resolution 
of the radio map is insufficient to correlate it with the X-ray 



blobs. However, there is some indication that, as observed 
in other clusters (see McNamara 2002 for a review), both X- 
ray holes may contain some steep spectrum radio emission 
which may be 'relic' emission from previous radio lobes. Fu- 
ture high resolution, low frequency radio observations will be 
crucial to linking the radio mini-halo and X-ray structure. 

It is worth noting that Fig.^]also suggests the presence in 
2A 0335+096 ghost X-Ray bubbles at larger radii similar the 
ones observed in other clusters (McNamara 2002). Unfortu- 
nately, because of the relatively short exposure, these bubbles 
are not statistically significant. A deeper X-ray observation is 
needed to confirm their existence. 

6. DISCUSSION 

The Chandra observation of 2A 0335+096 reveals that the 
core region of this otherwise relaxed cluster is very dynamic. 

The cluster image shows an edge-like feature at 40 — 90 kpc 
from the X-ray peak to the south. The edge spans a sector 
from 1 50° to 2 1 0° from the North. 

On smaller scales (< 20 kpc) the cluster image shows a 
very complex structure with X-ray filaments or blobs. More- 
over it shows the presence of two X-ray cavities (see Fig. [5}. 

The spectral analysis clearly shows that on large scales the 
clust er tem pera ture distribution is not azimuthally symmetric 
(see £14.21 and £|4.3I >. Furthermore, the cluster gas tempera- 
ture reaches its maximum in a region at « 90 kpc from the 
X-ray peak to the south, just beyond the surface brightness 
edge (Fig.|6jl. On smaller scales the X-ray blobs have similar 
temperatures but different pressures. 

Below we discuss the results of our analysis and propose 
possible dynamical models that may explain the observed 
properties of the core of 2A 0335+096. 

6.1. Nature of the surface brightness edge: a cold front 

In iQ.2l we analyze the spatial structure of the X-ray front. 
We show that the front is well fit by a discontinuous density 
profile with a density jump at the front of a factor i=s 1.6. Sim- 
ilar density jumps have been observed by Chandra in many 
clusters, including A2142, A3667, RX J1720+26, A1795. In 
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those clusters, they are interpreted as sharp boundaries of a 
dense, cooler gas cloud moving with a speed v through the 
hotter and more rarefied ambient gas. The surface bright- 
ness edge in 2A 0335+096 appears to be a similar cold front. 
As discussed by Vikhlinin, et al. 2002a, the pressure profile 
across the front can be used to determine the Mach number 
M = vi/cgi of the moving cooler gas cloud with respect to 
the speed of sound in the gas upstream, c s \. 

In Table 2 we show that the pressure jump across the front 
is Ap — 1.6 ± 0.3. This corresponds to a Mach number for 
the dense central gas cloud of M ~ 0.75 ± 0.2 (see e.g. Fig. 6 
of Vikhlinin, et al. 2001a). Such subsonic motion of the cool 
gas cores appears to be very common among the otherwise 
relaxed clusters with cooling flows (Markevitch, Vikhlinin, & 
Forman 2002). 

6.2. X-Ray Blobs 

In the analysis above, we identified a number of X-ray blobs 
most of which are likely to be part of a filamentary structure. 
As can been seen from Table 3, the blobs are quite dense: 
their densities range from 0.06-0.13 cm~ 3 . Although the 
projected spectral properties of the blobs are significantly dif- 
ferent (i.e. they have different emission weighted tempera- 
tures; see right panel of Fig. [6] and Fig.|9jl, in £14.51 we show 
that the temperature of their immediately surrounding ambi- 
ent gas and the deprojected blob temperatures are consistent 
within the errors for all but one of the blobs (see upper panel 
Fig. I10i. This is particularly interesting as it implies that the 
thermal pressure of the blobs is higher than the thermal pres- 
sure of the surrounding ambient gas. This finding, however, 
does not necessarily mean that the blobs are out of pressure 
equilibrium with the ambient gas as the latter may have extra 
non-thermal pressure (for example, magnetic fields and/or rel- 
ativistic electrons ejected from the central AGN). Although at 
the moment it is not clear which physical processes may be re- 
sponsible for the formation of the observed substructures, we 
can safely state that the system is dynamically unstable and 
the presence of blobs is likely to be a transient phenomenon. 
We notice that the evolution of the system strongly depends 
on the total pressure supporting the ambient gas. Below we 
discuss two extreme possibilities: 

1) the ambient gas surrounding the blobs is only supported 
by thermal pressure and thus the blobs are actually not in pres- 
sure equilibrium; 

2) the ambient gas, and possibly to some extent the blobs 
too, are supported by an extra non thermal pressure compo- 
nent such that the blobs and the ambient gas are actually in 
pressure equilibrium. 

In the first case, the blobs are overpressured so we expect 
them to expand. The expansion time is of the order of the 
blob's sound crossing time: 



tr 




(jj/0.69)m p 
(7/1.66)(/cT/1.3koV) 



(6) 



where r\ is the radius of the blob, and /j,, m p , and 7 are the 
mean particle weight, the proton mass and the adiabatic index, 
respectively. Using the linear dimension of the blobs reported 
in Table 3 we find t rj 1 — 2 x 10 7 yr, which appears to be very 
short compared with the likely cluster age. Furthermore, the 
blobs cooling time is w 10 times longer than the dynamical 
time (see lower panel of Fig. 1101 . Thus, in this first case in 
which the external gas is not supported by an extra pressure 



component, the blobs should expand before their excess pres- 
sure is radiated away. Consequently we do not expect to have 
a significant mass deposition into a cold phase from the blobs. 

Conversely, if the blobs are both in thermal and pressure 
equilibrium with the ambient gas, as in the second case con- 
sidered above, the blobs are likely to be more stable. Nev- 
ertheless, as the blobs are heavier than the ambient gas (they 
have higher densities), they are expected to sink eventually 
toward the minimum of the potential well. The time scale t s 
needed by the blobs to reach the cluster center is obviously 
t s >> tff 9 , where tff is the free fall time. If we assume the 
cluster potential well is described by a NFW profile (Navarro, 
Frenk & White, 1996) the free fall time for a point mass at 
rest at a radius r << r s can be well approximated by: 



(7) 



where r s , is the scale radius, g — 2irG5 c p cr u, G is the New- 
ton's constant, p crit — 3H 2 /8ttG is the critical density for 
closure, and 8 C is the characteristic density contrast. Using 
the best fit NFW parameters for our cluster and r w 20 kpc 
we find: 



t 



ff 



108 x hn 



20 kpc 



(8) 



From Table 3 we see that the free fall and the cooling times 
are similar, thus the actual sinking time is much longer than 
the cooling time. Therefore, unlike the first case considered 
above, if the ambient gas is supported by an extra non ther- 
mal pressure component such that the gas and the blobs are in 
pressure equilibrium, we expect the blobs to cool down before 
they can actually sink into cluster center. Hence, some mass 
deposition into a cold phase is possible. 

6.3. Possible dynamical models 

We conclude this section by discussing possible dynamical 
models that may explain the origin of the observed complex 
structure in the core of 2A 0335+096. Our starting point is 
the evidence for the presence of a dense gas core that moves 
from North to South, with a Mach number M ~ 0.75 ± 0.2 
(see 3Q . 

We propose two models: 

i) the motion of the cool core induces instabilities that pene- 
trate inside the core and disturb the gas, 

ii) the central galaxy has an intermittent AGN which makes 
the cool gas "bubbling". 

The features observed in the cluster result from the develop- 
ment of hydrodynamic instabilities induced by these two phe- 
nomena. Below we discuss in detail these two models. We 
conclude that both models are consistent with the available 
data and that further observations are required to discriminate 
between the two cases. 

6.3.1. Core Destruction by Hydrodynamic Instabilities 

We first show that small-scale structure can be the result of 
hydrodynamic instabilities induced by the observed motion 
of the core gas (revealed by the presence of the cold front. 
This motion may be due to a merger (as in A2142 or A3667, 
Markevitch et al. 2000; Vikhlinin et al. 2001), or due to the gas 
sloshing as in many other cooling flow clusters (Markevitch et 

9 Note that t a = tff just in the extreme case in which we assume that 
both the gas is collisionless and the blob have zero angular momentum. 
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al. 2001, 2002). The precise nature of the core motion is not 
important. 

The orientation of the front (see e.g. Figll 3i. together with 
the asymmetry in the temperature distribution (see FigHJ, 
suggests that the core is moving along a projected direction 
that goes from north to south. If the core is in fact a merging 
subcluster, it is clear that in projection the posit ion o f the sub- 
cluster is very close to the cluster center (see £13 . II . from the 
available data, however, we cannot measure the merger im- 
pact parameter. Hence, the subcluster may in fact be passing 
through the cluster at any distance along the line of sight from 
the center. 

When a dense gas cloud is moving with respect to a rarefied 
one, the interface between the two gasses is subject to hydro- 
dynamical Rayleigh-Taylor (R-T) and Kelvin-Helmholtz (K- 
H) instabilities. In the following we estimate the effect of 
both instabilities on the moving core using the same approach 
as Vikhlinin & Markevitch (2002). We assume that the gas 
cloud is a sphere of radius equal to the distance of the ob- 
served front, R ss 56.7" (« 39.7 kpc). As the cloud moves, 
the external gas flows around its border. The inflowing gas 
slows down at the leading edge of the sphere but reaccelerates 
at larger angles as it is squeezed to the sides by new portion 
of inflowing gas. In the case of an incompressible fluid, the 
velocity at the surface of the sphere is purely tangential and 
is given by v fi u id = 3/2^00 sin ip, where Voo is the veloc- 
ity of the inflowing gas at infinity and ip is the angle of the 
considered point on the cloud with respect to its direction of 
motion (see e.g. Lamb 1945). Unfortunately there is no an- 
alytic solution for the flow of a compressible fluid around a 
sphere. Nevertheless, the qualitative picture is similar. As we 
are only interested in order of magnitude estimates, in the fol- 
lowing we parametrize the actual fluid speed introducing the 
parameter rj so that: 



v fluid = Woo sin tp. 



(9) 



The parameter ?/ should range between 1 (fluid velocity at the 
surface of the sphere equal to the inflowing gas at infinity) and 
1.5 (incompressible fluid). 

First we consider the R-T instability. This instability de- 
velops only if the acceleration induced by the drag force d g 
is greater than the gravitational potential acceleration near the 
front due to the cluster mass inside the moving cloud. If we 
assume hydrostatic equilibrium, then the gravitational accel- 
eration is given by: 



9 = 



kT ( ci log n dlogT 



fiiripR \G?logr dlogr 



(10) 



where, {j,m p is the mean particle mass, and n and T are the 
gas electron density and temperature profiles, respe ctive l y. In 
our case this gives g = 1.7 x 10~ 8 cm s -2 (see £13 .21 ^4.4-1 
and Fig.|4j. On the other hand, the cloud acceleration due to 
the drag force is given 



g d = 0.15 



flout R 



(11) 



(where rii n ,n ou t the gas electron density inside and outside 
the edge) which in our case gives <jy = 2.5 x 10~ 9 cm s~ 2 . 
As the drag acceleration is less than the gravitational one, then 
the R-T instabilities should be suppressed. 

Now we consider the K-H instability. This instability may 
develop whenever two fluids in contact have a non-zero tan- 
gential velocity. As for the R-T instability, the presence of a 



gravitational field may suppress the development of the insta- 
bility. The stability condition is such that a perturbation with 
wavenumber less than k = 2ir/\ is stable if: 



Dk{v) 2 
D 2 - 1 



(12) 



where g is the gravitational acceleration at the interface be- 
tween the two regions and D = n in /n ou t is the density con- 
trast. Using the Eq.|9]for the tangential velocity, Eg. 1121 be- 
comes: 



A > 2?r- 



D r) sin ip 



400 - 900 sin 2 <p kpc. (13) 



D 2 - 1 g 

Eq. shows that gravity suppresses instabilities only on 
scales larger than 400 kpc. Hence, in absence of other sta- 
bilizing mechanism, we expect the formation of the K-H in- 
stabilities on smaller scales. 

In the case of 2A 0335+096 , we see the instabilities on 
scales of the order of the radius of the sphere R may de- 
velop already at ip > 18°. The condition that the instabil- 
ity on the gas cloud scale can develop, although necessary is 
not sufficient. In fact, in order to effectively disrupt the gas 
of the moving gas cloud, the growth time r of the instabil- 
ity must be much shorter than the dynamical time of the sys- 
tem (the cluster crossing time, t cross ). We thus require (see 
Vikhlinin, Markevitch, & Murray 2001b, and algebraic cor- 
rection of Mazzotta et al. 2002): 



tr. 



L 

3.3— simp >> 1. 
A 



(14) 



If we assume that the cluster radius is r = 1 Mpc, at the angle 
at which the instability starts to develop of ip = 18°, we find 
that t cross /T w 25. Furthermore, at the much larger angle of 
tp = 90°, where the instabilities are most effective, we find 
tcross/T ~ 83. Therefore, they have sufficient time to grow 
to the non-linear regime. 

The effect of K-H instabilities is to turbulently mix the gas 
in the dense moving cool cloud with the more diffuse hot- 
ter cluster gas. In particular some less dense gas from the 
cluster could be deposited into the denser gas cloud and vice- 
versa. Such a process may be responsible for the formation 
of the observed X-ray clumpy structure. Similarly, any pre- 
existing cold (< 30 K) gas deposited in the cluster and/or 
sub-cluster core by a cooling flow will also be mixed and re- 
distributed on these scales. This would explain the observed 
Ha+[Nll] emission morphology. On the other hand, these 
sub-sonic K-H instabilities will not produce any significant 
radio signature. Hence, in this scenario the radio properties 
of 2A 0335+096 (the mini-halo and the X-ray holes) must be 
related to activity of the central galaxy. 

In conclusion, the observed X-ray properties of the 2A 
0335+096 core may be produced by K-H instabilities induced 
by the motion of the core. 

6.3.2. Gas "bubbling" by central AGN 

The observed structure of 2A 0335+096 can also be ac- 
counted for if the central cluster galaxy hosts an AGN that 
undergoes intervals of strong activity followed by periods of 
relative quiescence (e.g. Binney & Tabor 1995, Soker et al. 
2000). During each active period, the AGN inflates two bub- 
bles by filling them with relativistic electrons. As observed in 
other clusters (e.g. Perseus, Fabian et al. 2000; Hydra A, Mc- 
Namara et al. 2000; Abell 2052, Blanton et al. 2001), this ex- 
pansion creates an external shell where the gas is denser. The 
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bubbles rise buoyantly, entrain the dense gas underneath them 
and induce convection (e.g., Churazov et al. 2001). These 
bubbles are also subjected to hydrodynamic instabilities that 
destroy them, fragmenting the denser shell into what we ob- 
serve as denser X-ray blobs. The observed structure in the 
core may be created by numerous bursts of AGN activity dur- 
ing which the jets point in different directions. We can vi- 
sualize this process as "gas bubbling" induced by the energy 
injection from a central AGN. The difference from other clus- 
ters with similar X-ray cavities but lacking such filamentary 
or blobby structure might be caused by a shorter duty cycle of 
the central AGN in 2A 0335+096, of the order of 10 7 yr (see 

The relativistic electrons produced by the AGN that are 
likely to be present inside the non-thermal bubbles (see e.g. 
McNamara 2002 and references therein) should eventually 
distribute in the central cluster region. However, they may 
avoid the densest gas blobs due to, for example, an enhanced 
magnetic field created during the shell compression induced 
by the bubble expansion. These electrons would add a non- 
thermal pressure component in the ambient (inter-blob) gas 
that may account for the apparent pressure nonequilibrium of 
the cluster core. 

As discussed in %.2\ in this case the dense X-ray blobs 
have time to cool, depositing some mass in a cold phase. 
This cool gas may explain the presence of the observed CO 
and Ha+[Nll] extended emission in the cluster center (see 
Fig.Hl. 

The relativistic electrons produced by the AGN may also 
account for the presence of the mini radio halo observed in 
the cluster core (see Fig. ll3l >. 

We finally mention that there may be a connection between 
the AGN activity and the core sloshing, hinted at by the ob- 
served correlation between the presence of the X-ray bub- 
bles and cold fronts in the central regions of many cooling 
flow clusters (Markevitch et al. 2002). Reynolds, Heinz & 
Begelman (2002) showed that during an AGN jet formation, 
a spherical pulse forms and propagates at the sound speed 
into the intracluster medium. Repeated AGN activity may in- 



ject sufficient kinetic energy to produce sloshing of the gas, 
although it is unclear how the bulk motion of the core as a 
whole could be produced. Sloshing of the dense core may 
also be induced by a merger, even if the merging subcluster 
did not reach the center of the main cluster (e.g., Churazov et 
al. 2003). 

To conclude this section, "gas bubbling" induced by central 
AGN provides a plausible explanation of the observed fea- 
tures in the cluster core. 

7. CONCLUSION 

The poor cluster 2A 0335+096 has several remarkable fea- 
tures that are only detectable with the superb resolution of 
Chandra. First, there is a cold front indicating a mildly sub- 
sonic gas motion. Second, two X-ray cavities that may be 
associated with steep spectrum radio emission indicating pre- 
vious radio activity in this system. Third, a number of small 
dense gas blobs in the cluster core that may be the shreds of 
a cooling core disturbed either by K-H instabilities or "bub- 
bling" induced by intermittent AGN activity. All of these 
properties relate to processes that may act to disrupt or destroy 
any cooling flow (sub-cluster merger, injection of non-thermal 
electrons into the intracluster medium or direct radiation from 
an AGN). From this relatively short X-ray observation it is not 
possible to determine which process (if any) dominates. Fu- 
ture deeper observations in the X-ray, radio and optical of 2A 
0335+096 may allow the energetics of each to be assessed. 
When combined with observations of other similar clusters, 
it will be possible draw concrete conclusions about the astro- 
physical nature of cooling flows. 
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